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Abstract 

This is the first part of lectures about effective field theories. Decoupling of 
heavy-particle loops is considered (heavy leptons in QED, heavy quarks in QCD). 
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We don't know all physics up to infinitely high energies (or down to infinitely small dis- 
tances). Therefore, all our theories are effective low-energy (or large-distance) theories 
(except The Theory of Everything, if such a thing exists). 

There is a high energy scale M (and a short distance scale 1/M) where an effective 
theory breaks down. We want to describe light particles (with masses mj <^ M) and their 
interactions at low energies, i.e., with characteristic momenta Pi <^ M (or, equivalently, at 
large distances ^ '^/M). To this end, we construct an effective Lagrangian containing the 
light fields. Physics at small distances < 1/M produces local interactions of these fields. 
The Lagrangian contains all possible operators (allowed by symmetries of our theory). 
Coefficients of operators of dimension n + 4 are proportional to If M is much larger 

than energies we are interested in, we can retain only renormalizable terms (dimension 4), 
and, perhaps, a few power corrections. 

These lectures can be used as an addendum to a standard textbook on quantum field 
theory. We follow notation of pQ ; details of some computations omitted here can be found 
in this book. Heisenberg-Euler Lagrangian is discussed in any QED textbook, see also 
the review [2]. Technically, the present lectures are largely devoted to QED. But the 
same ideas and methods are also used in QCD, only calculations are more lengthy: more 
diagrams, colour factors. . . The most clear treatment of decoupling by matching on-shell 
matrix elements in the full theory and the effective theory (where all loop corrections vanish 
because they contain no scale) is presented in [3] . The method of regions, which is a useful 
alternative to constructing effective Lagrangians, is discussed in the textbook 

2 Heisenberg— Euler effective theory 
2.1 Photonia 

The first historical example of an effective low-energy theory is the Heisenberg-Euler ef- 
fective theory in QED. It is still the best example illustrating typical features of such 
theories. 
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In order to understand it better, let's imagine a country, Photonia, in which physicists 
have high-intensity sources and excellent detectors of low-energy photons, but they don't 
have electrons and don't know that such a particle exist£]. At first their experiments 
(Fig. [1^) show that photons do not interact with each other. They construct a theory. 
Quantum PhotoDynamics (QPD), with the Lagrangian 



(2.1] 



But later, after they increased the luminosity (and energy) of their "photon colliders" and 
the sensitivity of their detectors, they discover that photons do scatter, though with a very 
small cross-section (Fig. [T]d). They need to add some interaction terms to this Lagrangian. 





Figure 1: Scattering of low-energy photons 
There are no dimension 6 gauge-invariant operators, because 



F/F/F/ = 



(2.2) 



(this algebraic fact reflects C-parity conservation). The operator F^^d'^F^^^ can be rewritten 



as 



F^.d'^F^^ = -2F^M''F''' = 2d^ {F^.dxF^") - 2 (d^F,^) {OxF-"^) , (2.3) 
where the identity 

dxF^, + d^F^x + d.Fx^, = (2.4) 

has been used. The full derivative can be omitted from the Lagrangian. So, only one 
dimension 6 operator with two F^^^ remains, 

O = (d^Fx,) {d.F^") . (2.5) 

But this operator vanishes due to equations of motion 



a.F 



0. 



(2.6) 



^We indignantly refuse to discuss the question "What the experimentahsts and their apparata are made 
of?" as irrelevant. 



On-shell matrix elements of such operators vanish; therefore, we may omit them from the 
Lagrangian without affecting the S'-matrix. 

Interaction operators first appear at dimension 8: 

Oi = {F,,Fn' , O2 = F.^F-'^^Fo^pF"^ . (2.7) 

Hence the QPD Lagrangian which incorporated the photon-photon interaction is 

L = Lo + Li, Li = ciOi + C2O2 , (2.8) 

where the coefficients ci^2 ~ l/M'^, M is some large mass (the scale of new physics). Of 
course, operators of dimensions > 8 can be also included, multiplied by higher powers of 
1/M, but their effect at low energies is much smaller. Physicists from Photonia can extract 
the two parameters Ci^2 from two experimental results, and predict results of infinitely many 
measurements. 

We are working at the order therefore, the photon-photon interaction vertex 

can appear in any diagram at most once. The only photon self-energy diagram vanishes: 

^^^^=0, (2.9) 

because the massless loop is scale-free. Therefore, the full photon propagator is equal to the 
free one. There are no loop corrections to the 4-photon vertex, and hence the interaction 
operators (12.71) don't renormalize. 

2.2 Qedland 

In the neighboring country Qedland physicists are more advanced. In addition to photons, 
they know electrons and positrons, and investigate their interactions at energies ~ M 
(M is the electron mass). They have constructed a nice theory, QED, which describe their 
experimental result^. 

Physicists from Qedland understand that QPD constructed in Photonia is just a low- 
energy approximation to QED. The coefficients Ci^2 can be calculated by matching: we 
calculate the amplitude of photon-photon scattering at low energies in both the full theory 
(QED) and in the effective theory (QPD), and equate them. In QPD the scattering am- 
plitude is T^^'^^-^-^(pi,p2,K,p'2)ei^,e2;.,e;;^e'2;^ where T'^i/^^-i-^ = c^Tj^i'^^'^^^^ + c^T^^'^"''''', 
and the tensors j^i^^^^^^'^ can be trivially obtained from the operators Oi^2 (12. 7p (they are 
linear in each momentum). 

In order to find 2 coefficients Ci_2, it is sufficient to match just two scattering amplitudes. 
For example, we can consider forward scattering: p'^ = pi, p'2 = P2, Pi = = 0, {pi + 

^They don't know muons, but this is another story (Sect. [3]). 
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= s <^ M^, where M is the electron mass. It is sufficient to calculate 2 numbers in 
QED: T^^^^'^^^"'^ g ^^y^g ^.^i,^ and T^^'^^'-'^^'^g^^f^^g^^^^, and then Ci^2 can be obtained by solving 
the linear system. In QED, photon-photon scattering ffist appears at one loop and is 
given by 6 diagrams; pairs of diagrams which differ from each other only by the direction 
of the electron line contribute equally (just insert Uq = 1 between all propagators and 
vertices, where Uc is the charge conjugation matrix: Ucl^Uc = — 7J, where T means 
transposition). Hence the scattering amplitude T'^^^^'^^'^^ {pi, p2, pi, P2) is twice the sum of 
3 diagrams (Fig. [2]). 
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Figure 2: Photon-photon scattering in QED at one loop 

In order to calculate the scalar quantities T'^^^^'^^^^'^ g ^^^^g ^^^^ and T^^'^'^^^^'^ g ^-^^^g^^^^, we 
expand the propagators in pi 2 <^ M up to the fourth order, and average the integrand in 
k directions. The results are expressed via the one-loop vacuum integrals (Fig. [3D 

-Lj I ^ = M--t/(") . D^NP-k^-m, V(n) = ^J^. (2.10) 
Solving the linear system for ci_2, we arrive at the scattering amplitude 

4 /If -4-2e /J _ A\(A _ a\ 

(47r)'^/2 ^ ^ 2880 ^ ^ ^2 ) ■ K ■ ) 




Figure 3: One-loop massive vacuum diagram 



5 



It is finite at e — > 0. To reproduce this QED result, the interaction term in the QPD 
Lagrangian fl2.8p must be [5] 



-50i + UO2) 



(2.12) 



180M4 

This effective Lagrangian can be apphed to any problems involving soft photons and 
their interactions. As an example, let's consider the thermal radiation at a temperature 
T <C M. What is its energy density? The leading term is given by the first diagram in 
Fig. m Note that this diagram is not scale-free and does not vanish, because T sets the 
scale. The energy density is ~ T^, by dimensionality (Stefan-Boltzmann law). The first 
correction is given by the second diagram in Fig. HI It contains the coupling ~ a^/M^, 
and, by dimensionality, it gives [6] 



a 



rp8 



It would be difficult to guess this result without using the effective theory. 





Figure 4: Energy density of thermal radiation 



It is not difficult to calculate the two-loop correction to the Lagrangian (12.121) . The 
two-loop scattering amplitude in QED reduces to the two-loop vacuum integrals (Fig. [5]) 



(m^/2)2 J DfwflDf 



Di = M^-kl, D2 = M^-kl, D^ = -{ki-k2) 



(2.13) 



V{ni,n2,n^) 



F (f - ng) F (m + ns - f ) F (^2 + ^^3 - f) r(ni +n2 + 71^-8) 



F (f ) F(ni)F(n2)F(ni + 713 + 2n^ - d) 
-iO is assumed in all denominators). 




Figure 5: Two-loop massive vacuum diagram 
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2.3 Coulomb potential 

Let's suppose that physicists in Photonia have some classical (infinitely heavy) charged 
particles, and can manipulate them at their will. If a particle with charge e moves along a 
world line /, the action contains the interaction term 



■'int 



(2.14) 



in addition to the photon field action. The integrand exp(z5') in the Feynman path integral 
contains a phase factor 



Wi = exp ( ie 



dx^ i^x ^ 



(2.15) 



called the Wilson line (when the operator language is used, the exponent is path-ordered). 
The vacuum-to-vacuum transition amplitude in the presence of classical charges is thus 
the vacuum average of the corresponding Wilson lines §. 

Suppose two charges e and — e stay at rest at some distance r during some (large) time 
T. The energy of this system is f/(r) — the interaction potential of the charges. The 
vacuum transition amplitude is 

e-^^(^^)^. (2.16) 

On the other hand, this amplitude is the vacuum average of the Wilson loop (Fig. [6]) with 
T ^ r (we don't care what happens near its lower and upper ends as long as these regions 
are small as compared to T). 



Figure 6: Wilson loop 

The zeroth-order term in the vacuum average of any Wilson loop is 1. It is convenient 
(though not necessary) to use the Coulomb gauge to calculate the first correction. In this 
gauge, there is Coulomb photon with the propagator 



q- 



(2.17) 



•^Later we'll see that the propagator of a heavy charged particle in the effective theory which describes 
its interaction with soft photons in the Wilson line. 
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(it propagates instantaneously) and transverse photon with the propagator 



1 



g2 + iO 



5'^ 



(2.18) 



Wilson lines along the direction only interact with Coulomb photons. The self-energy of 
the classical particle vanishes: 



0, 



because the particle propagates along time, and the Coulomb photon along space. 
Therefore, there is just one contribution at the order e^: 

T 



(2.19) 



r 



,',r + 1 







le'T I D^\t,r)dt = -ie'T 



" y noo 



(2vr; 



(2.20) 



(integration in r gives T). Comparing it with 1 — iU{f)T (12.161) . we obtain the Fourier 
transform of the potential 

[/(g)=e^D°°(0,g) = -4; 
at (i = 4 the Coulomb potential is 



f2.211 



U{t\ 



a 
r 



(2.22) 



What about corrections? Vertex corrections (Fig. [Th) vanish for the same reason 
as (12.191) . Crossed-box diagrams like in Fig. [7^ vanish because Coulomb photons prop- 
agate instantaneously, and the time orderings of the vertices on the left line and on the 
right one cannot be opposite. We don't need two-particle-reducible diagrams like in Fig. ^ 
because they match higher orders of expansion of the exponent (12.161) . Only corrections to 
the photon propagator can contribute. But there are no such corrections in QPD. Hence 
the Coulomb potential (12.221) is exact in this theory. 

In the presence of sources of the photon field, the dimension 6 operator O (12. 5p cannot 
be ignored. The QPD Lagrangian now contains an extra term 



cO 



(2.23) 
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X 



Figure 7: Corrections to the Wilson loop 

where c ~ by dimensionality. This term produces the contribution to the photon 

self-energy 

'^v^wwBvw^ = 2icq^ {q^g,,v - Qf^qu) ■ (2.24) 
q q 

The aim of an effective theory is to reproduce the S-matrix, so, we may neglect operators 
vanishing due to equations of motion (in particular, the self-energy (12.241) vanishes at the 
photon mass shell = 0, but is needed in a virtual photon line exchanged between classical 
sources). Therefore, the term (12.231) reduces to cj^j^, where is the external (classical) 
current. For the classical charges e and — e, it leads to the contact interactioij^ 



U,{r) = 2c5{r) . 



(2.25) 



What can physicists in Qedland say about the interaction potential between classical 
charged particles? The photon self-energy in QED at one loop is 



k + q 




[q^dtiv - qt^qv) n(g^ 



(2.26) 



It is convenient to contract this with g^^: 



n(g^ 



Aiel 



{d-l)q^J {27rYD,D2' 



Ml 



{k + qf 



Do = M, 







iV = - Tr 7^(1^ + ^ + Mo)r{lt + Mo) 



'-(Di + D2 + q' 



2Ml 



(2.27) 



(here Cq is the bare electron charge, and Mq is the bare electron mass). When q -C M, we 
can expand 1/Di in g; then averaging over k directions reduces the problem to the scalar 



"'it is difRcult to observe a i5-function potential in the interaction of classical charged particles. But this 
interaction is essential if the particles are quantum-mechanical — it shifts energies of S'-wave states (we'll 
discuss this later, in NRQED). 
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integrals fl2.10p . The result is 



-2e 



3 (47r)'^/2 



T(e) 1 



10 M2 



+ 



f2.28) 



If we retain only the leading term in fl2.28p . then the potential U{q ) in QED has the 
Coulomb form fl2.2ip . From now on, the charge in the low-energy effective theory (QPD) 
will be denoted e' (it is not renormalized: e' = Cq). The interaction potentials in both 
theories match if 



J2 



AeiM: 



-2e 



w/o r(£) . (2.29) 
3 (47r)'^/2 ^ ^ ' 

The renormalized QED charge e(/x) in the MS renormalization scheme is related to the 
bare one Cq by 



e (/i)Zc,(a(/i)) 



a 



Z«(a) = l-/3o— + 



(2.30) 
(2.31) 

47r [ATcy/^ ^ ^^-^^^ 

(7^; is the Euler constant). Let's express the running QED charge a{fi) via the QPD charge 
e' from fl2.29p with the one-loop accuracy: 



where 



and 



a(/i) e2(yu) At 



-7-B£ 



47r 



(47r)^/2 



1 + 



-2e 



/5o /M\ 



2e 



(47r)^/2 \ e V / 



— + -e^^T(£) + ■ 



(2.33) 



In the correction term we may replace Mq by a renormalized electron mass M (in any 
scheme — differences are beyond our current accuracy). Both e' and a{fi) are finite at 
£ — > 0. Therefore, 

Po = -t- (2-34) 



The running QED coupling a(/i) in the physical limit e — > is expressed via the QPD 
a' = e'2/(47r) (which is measured in macroscopic experiments) by 



a a 
a(/i)=a'(l-/?o— log— + 



It satisfies the renormalization group equation 

d loga(/i) 



a 



d log fi 

(see f l2.34p ) with the initial condition 



-2/3(a(/i)), (3{a) = l3o— + 

Air 



(2.35) 



(2.36) 



a{M) = a' . 



(2.37) 
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What about the term in (12.281) ? Comparing it with fl2.24p we obtain 



at e = 0, and thus the contact interaction fl2.25p . We can also see this more directly: this 
term cancels the denominator of the photon propagator in (I2.2ip . and thus give the 
contribution 

to the interaction potential at e = 0. 

2.4 The full theory (QED) and the low-energy effective theory 
(QPD) 

Now we shall discuss relations of these two theories more systematically. The Lagrangian 
of the low-energy effective theory (QPD) is 

^' = -\Po,uFr - ^^{d,Kf (2.40) 

(primed quantities are those of QPD; we neglect the quartic interaction terms for now). 
There are no loops which can be inserted into the photon propagator in this theory, and 
hence no renormalization: 

A'o = A'(/i) = ^',3 , a'o = a'(/x) = (2.41) 

(the index "os" means the on-shell renormalization scheme; quantities with the argument 
11 are those in the MS scheme). 

The QED bare fields and parameters are related to renormalized ones in the MS scheme 

by 

V^o = ^y^(a(/i))V^(/i) , Aq = Z^J\a{fi))A{fi) , 

ao = ZA(a(Ai))a(/i) , Cq = Zy2(a(/i))e(/i) , Mq = Zjn{a{^j))M{^) , 
where all renormalization constants have the minimal form 

-7i^ + (f + f)(a^- 

and a(/i) is defined by (12.321) . 

Another renormalization scheme often used in QED is the on-shell scheme. The photon 
field renormalized in this scheme is related to the bare one by 

Ao = (Z°^)'/'Ao, (2.44) 
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(where the renormahzation constant Z'^ is not minimal), and therefore the bare and renor- 
mahzed photon propagators are related by 

D^ip') = Z^^DTip') . (2.45) 
The photon propagator near the mass shell is 

By definition, Df{p'^) behaves as the free propagator l/p"^ near the mass shell. Therefore, 

Zf = ^— . (2.47) 

^ l-n(O) ^ ^ 

We want to calculate processes with low-energy photons in QED, and compare its 
predictions with those of the low-energy theory, QPD. First of all, we shall consider the 
photon propagator. The propagators of both A^s and A'^^ at ^ are equal to the free 
propagator, and therefore 

Aos = Ks- (2.48) 
The bare photon fields in the two theories are related by 

Ao = iClY^" A', (2.49) 

(up to corrections suppressed by powers of 1/M). From (I2.48P we obtain the bare decou- 
pling coefficient 

C^ = z^^ (2.50) 

The corresponding relation between the MS renormalized fields is 

Aifi) = efif^)A'if^) , (2.51) 
where the renormalized decoupling coefficient is 

a(/i) = ^ = ^ • (2.52) 



At one loop 



(c^)" = {^ir' = 1 - n(o) = 1 + f^f^m + ■■■ (2.53) 



(see ( I2.28p . We can re-express it via renormalized quantities: 

|£^r,.)^e-e^T(l + .)^Z„Z-, i^21og^. (2.54) 
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The renormalized decoupling coefficient = Za{(a)^^ — ^a/Z^ niust be finite at £ 
(for example, at L = 0). We substitute the MS renormalization constant Za = 1 
Zia{fi)/{ATTe) with an unknown zi, and find from this requirement 



Finally, we arrive at 



G'W = 1 + + • ■ • (2.56) 

Calculation of n(0) at two loops (Fig. [8]) reduces to the vacuum integrals (12.131) . and 
it is not difficult to obtain the result 

n(0) - -^'l^r(e) - 2 (^-4) (5^^ -33^ + 34) (ejM^^.y ... 





Figure 8: Two- loop photon self-energy 

Expressing the on-shell renormalization constant (12.471) (and the bare decoupling coef- 
ficient (I2.50p ) via renormalized quantities (I2.54p we obtain 

(CS)"' = (z?)-' = 1 - n(o) = 1 + li=^z„z-^- 

n ^ ^ '2,58) 



When re-expressing it via renormalized quantities, we have to use one-loop renormalization 
constants in the one-loop term. The charge renormalization constant can be obtained from 
the photon-field one (12.551) : 

Z^ = Za' = ^ + 1^ + --- (2.59) 

The MS mass renormalization constant is 

Z^ = l-3^ + --- (2.60) 

(see (IA.17P ). The renormalized decoupling coefficient (^^ = Za/Z'^ must be finite at 
e ^ (for example, at L = 0). We substitute the MS renormalization constant Za = 
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1 — (4/3)a;(/i)/(47r£) + (^20 + Z2ie){a{fi) / {4:TTe)y with unknown coefficients Z20, Z21, and 
find from tliis requirement 

4 Q;(/i) / 

Z;^ = i - 

Finally, we arrive at 



Za=1 — - 2e +■■■ (2.61) 



Now we are going to compare the charge in QED and in the low-energy theory (QPD). 
In QED, due to Ward identities, 

Za = Z^' , = {Z^r' ■ (2.63) 

The on-shell charge (measured at large distances) is the same in both theories: 

aos = • (2.64) 

Therefore, 

= CMa\ii) , (2.65) 

where 

CM = ^ = Ca\i^)- (2.66) 

At what scale should we do this matching? In principle, this does not matter. 
However, in practice, we know only a few terms in perturbative series. Their coefficients 
contain powers of the logarithm L f l2.54p . If it is large, truncating the series after some 
term produces large errors. Therefore, it is better to use /i ~ M. One popular choice is /iq 
defined as the root of the equation 

yuo = M(/io) (2.67) 

(it corresponds to L = 0). Then (12.621) 



13 / a{iiQ) 



C(/Xo) = 1 + y (^^J +■■• (2.68) 

Another popular choice is to use the on-shell mass: /i = M^s- It is related to the MS mass 
by 

^^^)-l-6riog^ + ?U + -- (2.69) 



Mos V 3; 47r 

and we obtain 

C(MJ = 1 + I5(^^^y + --- (2.70) 

In general, for any /i which differs from the mass M by 0(a), Ca{^^) differs from 1 by 0{a'^) 
(this is so, in particular, for /io and Mqs). If we choose, for example, = 2M or = M/2 
instead, we'll have C,a{^^) = 1 + 0{a), i.e., the correction will be much more important. 
Therefore, it is better to match the full theory and the low-energy one at some scale 
which differs from Mqs (or /^o) by 0{q). 
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3 QED with heavy muons 



3.1 Qedland 

Physicists in Qedland suspect that QED is also only a low-energy effective theory. We 
know that they are right, and muons exislEI. There are two ways in which they can search 
for new physics: 

• by increasing the energy of their e^e~ colliders in the hope to produce pairs of new 
particles; 

• by performing high-precision experiments at low energies (e.g., by measuring the 
electron magnetic moment). 

New physics can produce new local interactions of photons, electrons, and positrons, which 
should be included in the effective QED Lagrangiaii^. 

For simplicity, we shall first assume that electrons are massless. This is a valid ap- 
proximation if we want to describe phenomena at energies much higher than nte but much 
lower than the new physics scale M (in reality, the muon mass). A non-zero electron mass 
will be re- instated in Sect. 13.81 What operators of dimensions > 4 can we add to the QED 
Lagrangian? The first gauge-invariant operator appears at dimension 5 (the anomalous 
magnetic moment): 

= ^F^,a^''^, (3.1) 

where a^^ = (i/2)[7'^, 7^]. However, it flips the electron helicity. Helicity of a massless 
electron is conserved in QED. If we assume that new charged particles interact with pho- 
tons, but not with electrons (or, if they interact with electrons, this interaction does not 
violate helicity conservation), this operator cannot appear in the Lagrangian at me = 0. 

Under this assumption, new interactions first appear at dimension 6. Electron-electron 
contact interactions 

On = (V'7(ri)V')(V^7(n)V^) (3.2) 

(see IA.2P conserve helicity at odd n. The operators 



reduce to Oi (13. 2p due to equations of motion, and 



^In our real world M.^ ~ M^; for simplicity we shall assume that pious don't exist. 

^We were lucky that the scale of new physics in QED is far away from the electron mass m. Contribu- 
tions of heavy-particle loops are also strongly suppressed by powers of a. Therefore, the prediction for the 
electron magnetic moment from the pure QED Lagrangian (without nonrenormalizable corrections) is in 
good agreement with experiment. After this spectacular success of the simplest Dirac equation (without 
the Pauli term) for electrons, physicists expected that the same holds for the proton, and its magnetic 
moment is e/(2mp). No luck here. This shows that the picture of the proton as a point-like structureless 
particle is a poor approximation at the energy scale m^. 
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Therefore, our first task is to investigate renormalization of the operators (13.21) . 

An important aspect of all effective field theories is power counting. In this particular 
case it is trivial; however, in more complicated situations, it becomes more involved. When 
considering QED processes with small characteristic momenta p, we have a small parameter 
A ~ p/M. The characteristic distance at which soft fields vary is large: a; ~ 1/A, and 9^ ~ A 
when acting on soft fields. The soft photon propagator is 



<0\T{A^{x)A,{0)}\0> j 



vp-x 

(27r)4 p2 



g^u - (1 - 



p2 



from p ~ A and a; ~ 1/A we obtain A ~ A (this means that ~ A is homogeneous in A). 
The soft electron propagator is 



<0\T {^{x)^{0)} |0> 



d P -ip-x ^ 



and we obtain if) ~ A^^^. The terms in the leading-order Lagrangian scale as F^^F^'^ ~ A^, 
i/jUpip ~ A^; this means that the characteristic action is of order 1, as expected. The first 
power corrections (13.21) to the Lagrangian scale as A^, and their contribution to the action 
is ~ A^. 

Of course, we can add higher-dimensional contributions to the Lagrangian, with further 
unknown coefficients. To any finite order in 1/M, the number of such coefficients is finite, 
and the theory has predictive power. For example, if we want to work at the order 
then either a single (dimension 8) vertex or two 1/M^ ones (dimension 6) can occur 

in a diagram. UV divergences which appear in diagrams with two dimension 6 vertices are 
compensated by dimension 8 counterterms. So, the theory can be renormalized. The usual 
arguments about nonrenormalizability are based on considering diagrams with arbitrarily 
many vertices of nonrenormalizable interactions (operators of dimensions > 4); this leads 
to infinitely many free parameters in the theory. As stressed already, in any effective theory 
we always work up to some finite order in 1/M, and the number of parameters is finite. 

3.2 Renormalization of four-fermion operators 

First we'll consider a somewhat simpler problem [7]: renormalization of the operators 

0° = (^107(n)^20)(V^307(n)^40) (3.3) 

with 4 different lepton flavours (7(„) is defined in (1A.2I) ). In 4 dimensions, only operators 
with n < 4 exist. However, in dimensional regularization bare operators with all n exist 
(those with n > 4 are called evanescent operators). If we use the standard minimal 
subtraction, renormalized evanescent operators don't vanish. Naturally, we want them 
to vanish, and thus to have only a finite number of renormalized operators. This can be 
achieved by tuning the renormalization constants (which are no longer minimal in the sense 

of dm). 
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a 



b 



c 






Figure 9: Matrix element of 



Let's calculate the matrix element of the bare operator (13.31) at one loop, keeping 
only the leading terms in the e expansion. We'll write the tree contribution (Fig. Maj) as 
7(n) ®7(n); it should be multiplied by the external legs renormalization factor (zV^)"^ (]A.9p . 



After averaging over k directions and using (lA.7p . the contributions of Fig. [9^, c; d, e; /, 
g are 



2 

-2 
+ 2 



■^l^lul{n)YY ® 7(n) - (1 - ahin) ® 7(n) 
^1^^7ulin) ® 7^7''7(n) - (1 " a)7(n) ® 7(n) 
^7(n)7!.7M ® 7''7''7(n) - (1 - a)7(n) ® 7(n) 



a 
47r£ 

a 
47re 

a 



The first of them is the same as for two-fermion currents, and has the structure 7(„,) ® 7(„) , 
see fraj) . 

The other ones can be calculated using the identities 



7^^7in) ® 7''7(n) = 7(n+l) ® 7(n+l) + ^(c? - n + l)7(n.-l) ® 7(n--l) , 
7(n)7M ® 7''7(n) = (-1)" [7("+l) ® 'yC^+l) - '^(C? - + l)7(n-l) ® 7(n-l)] 



(3.4) 



'''Here the zigzag line is actually a point; it only shows which fermion legs are connected to the same 7 
matrix. 
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Indeed, 

7(n) ® 7(n) = nl (7o7i ■ ■ ■ 7n-i ® 7°7^ ■ ■ ■ 7""^ H ) 

(there are d\/{n\ [d — n)\) terms in the sum). Now let's muhiply the left part (before ®) 
by 7^ and the right one by 7^^, say, from the left. There are two kinds of contributions. 
When the value of the index /i is not one of those already present, we get n + 1 7 matrices 
both before and after ®. Each term is produced n + 1 times; this factor converts n! into 
(n + 1)!, and the sum of all such contributions is 7(„+i) ® 7(n+i), the first term in the 
right-hand side of the first identity. Otherwise, one 7 matrix gets squared, and we get 
n — 1 7 matrices both before and after 0. There are n such values of fi for each term; the 
resulting bracket contains d\/{{n — l)\ (d — ra + l)!) terms, so, each one is produced d — n + 1 
times. Converting the common factor n\ to {n — 1)!, we get the second term in the first 
identity. The second identity follows easily: contributions of the first kind get an extra 
factor (—1)" from commuting 7^ to the left; contributions of the second kind get (— 
Note a very important property: when n = 5 (Dirac structure vanishing in 4 dimensions), 
the contribution 7(ra_i) <S>7(n-i) (Dirac structure not vanishing in 4 dimensions) comes with 
the factor ~ e. 

Combining all contributions, we obtain a gauge-invariant result for the matrix element 
of the bare operator (13. 3p : 



<0°> 



1 + 2(n- l)(n-3) 



a 
Aire 



7(n) ® 7(n) 



7(n+2) ® 7(n+2) + n{n - l){d - u + l){d - u + 2)7(„_2) ® 7(„-2) 



a 



(3.5) 



J 4:7ie 



In particular. 



/ <0?> \ 



<03°> 

<ol> 



\ 



J 



1 


-1 




\ 






/ 


7(1)^ 


5 7(1) 


\ 


-36 





-1 




a 






7(3) <^ 


5 7(3) 






40£ 


16 


-1 


4:7ie 






7(5)^ 


5 7(5) 




v 






■•• ) 






v 






/ 



(3.6) 



(where the physical operators with n < 4 and the evanescent ones (n > 5) are separated 
by lines). Here only the leading term in e is given for each element of the matrix. These 
leading terms originate from the UV divergent loop integral (]A.7|) : they don't depend on 
external momenta and masses. The contribution 7(3) ®7(3) in the matrix element of the bare 
evanescent operator O5 acquires the factor e from (13.41) . and hence is finite, not 0{l/e). 

The column vectors of the bare and renormalized operators are related by the matrix 
of renormalization constants: 



0° = Z(a(/i))0(/i) , O(^) = Z-\a{n))0^ . 
The renormalized operators obey the renormalization group equations 

dO{n) 



dlogn 



+ 7(a(yu))O(/i) = 0, 



(3.7) 



(3. 
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where the anomalous dimension matrix is 

_i dZ 

7 = 



dZ- 



-Z . 



d log /i d log /i 
We want the renormalized evanescent operators to vanish: 

0(^) = 



(3.9) 







V 



(3.10) 



Therefore, we have to include this 0{1) term in Z: 

( 

Z(a) = 1 + 



-1 




\ 




-36 


-1 




a 


40e 


16 


-1 
■•■ ) 


47re 



(3.11) 



Then 



05(/i) 



01 



and the 7(3) ® 7(3) contributions of 0° and 0° cancel in matrix elements of O^^fi) = 0. 
The renormalization constant matrix ( 13.1ip is not minimal: it contains not only negative 
powers of e, but also C(l) contributions. The one-loop anomalous dimension matrix is 



7(a) 



/ 


-1 


\ 




-36 





-1 


a 






16 -1 


4:71 


v 




'■■ ) 





(3.12) 



Its lower left block vanishes; therefore, the form fl3.10p is preserved by evolution (13.81) . 
Evolution of the physical operators is determined by the upper left block. When calculating 
this block at one loop, one may forget about evanescent operators. 

With the two-loop accuracy, the non-minimal renormalization matrix can be written 

as 



Z(a) 



no 



a 



An 



21 



Z 



22 



— ^+U20 + — + -f T- 



47r/ 



(3.13) 



The anomalous dimension matrix (13.91) must be finite at e 
consistency condition 

^22 = 7:^11(2^11 — Po) , 



0. This gives the self- 
(3.14) 



i.e., the l/e"^ part of the two-loop contribution Z22 is determined by the l/e part of the 
one-loop contribution Zu. The anomalous dimension is 

2 



7(a) 



a 



'2^11- 2(2^21 

47r 



/ a 

ZioZii — ZiiZiQ + PoZio) y— 



(3.15) 
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As we have already discussed, the one-loop 0{e ) term has the structure 

Z.o=(°°); (3.16) 

it originates from the product ofl/e divergences of one-loop integrals (which are momentum- 
independent) and the factor e from 7-matrix algebra. Because of this factor, there are no 
l/e terms in the lower left corner: 

Zn=(t (3.17) 



The lower left corner of 

Z21 



d 

e f 
g h 



comes from the product of divergences of two-loop integrals (which are momentum- 
independent) and the factor e from 7-matrix algebra. These l/e"^ divergences of two-loop 
integrals are determined by products of l/e divergences of one-loop integrals. Namely, 
the lower left corner of the anomalous dimension matrix fl3.15p must vanish, in order to 
preserve the form (13.101) . This leads to the self-consistency condition 

g = -{ab + da — Poa) . (3.18) 

Evolution (13. 8p of the physical (non-evanescent) operators is determined by the upper left 
corner of the anomalous dimension matrix 7 (13.151) : the two- loop contribution to it is 



-2 2e + ca — 
V47r 



a y 



In order to find it, one needs e, the l/e part of two-loop diagrams with the insertion of 
a physical operator. Such a calculation (though in a more complicated case) is explained 
in [8]. The best way is to nullify all external momenta, and to insert masses into all 
denominators to regularize IR divergences. 

In the case of the operators (13.21) with identical flavours, an electron line leaving the 
operator vertex can returns to the same vertex. The e"'"e~7 matrix element at one loop is 
given by the diagrams in Fig. [TUJ These diagrams contain the electron loop 

which vanishes after contraction with p^. Therefore, this matrix element is proportional 
to the tree-level matrix element of the operator O = ipd^F^^^'^^il). The first diagram 
contributes only to <0°>. Taking the trace of (13.191) with a single 7 matrix (and the 
fermion loop factor —1) we recover the familiar expression for the one- loop photon self 
energy. The divergent part of the matrix element is 
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Figure 10: Matrix element of 0\ 



(we have taken into account that there are 2 such contributions). 
The second diagram for <0°> contains 

7(n)7''7{n) = ^nl^ ■ 



(3.21) 



Multiplying this by 7^ and using flA.6p we have 



-ir 



d-ln 



d 



-l{n)l(n) 



There are d\/{n\{d — n)!) ways to select n indices of d, and they contribute equally to 

7(n)7(n): 

7(n)7(.) = ;^y^^7o---7n-i ■ ■ ■ 7""^ ± " • • ) ■ 

There are n! terms in the bracket; they contribute equally, namely, (— l)"("^i)/2 each. 
Finally, we obtain 



An = (-l)"("+^)/2(^_2^)^ 



T{d-n + l) 

(_^y«(n+i)/2(^ _2n).{d-l){d-2)---{d-n + l) 



(3.22) 



(naturally. An contains {d — 4) when n > 5). Using this algebraic fact, we arrive at the 
divergent part of the two contributions given by the second diagram in Fig. (TO) 



.(_l)«("+l)/2^ 



n-2 



{A-n)\ (47r)% 



<0>. 



(3.23) 



Due to the equation of motion, O = eOi. Collecting all contributions together, we 
obtain 



<0^,> 

<oi> 



1 + 



-4 -1 \ a 
-37 J 4^ 



<0»> 

tree 

^'-^3 ^tree 



21 



where only 1/e parts of the one-loop matrix elements are retained. This gives the anomalous 
dimension matrix 

7 = 2^ ,t l)^ (3.24) 

of the operators (13. 2p . 
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3.3 Contact interaction of electrons 

The Lagrangian of the effective low-energy QED (without muons) with the accuracy 
is 

L = Lo + Li , Li = c?0° + c°0° = ci(/i)Oi(/x) + Cs{fx)03{fx) , (3.25) 
where Lq is the standard massless QED Lagrangian, the bare operators 0° are defined by 

On = (V^07(n)V^o)(V^07(n)^o) , (3.26) 

and the coefficients ~ play the role of charges at the vertices produced by these 

operators (similarly to Cq at the ordinary electron-photon vertices). 

Interaction terms (of some dimensionality) in any effective Lagrangian can be written 

as 

Li = = c^(/i)0(/i) , Oo = ^(a(/i))0(/i) , c(/i) = Z^(a(/i))co , (3.27) 

where Oq is the column vector of the bare operators, Cq is the row vector of the bare 
couplings (T means transposition), and Z is the matrix of renormalization constants of 
these operators. The renormalized couplings (called also Wilson coefficients) obey the 
renormalization group equations 



(ic(/i) y 



dlog/i 



Y («(/i))c(/x) , (3.28) 



where the anomalous dimension matrix of the operators O is defined by (13. 9p . The Wilson 
coefficients Ci{^o) are determined by matching — equating some S'-matrix elements in 
the full theory (expanded in Pi/M) and in the effective theory. It is most convenient to 
use some /iq ~ M; then Cj(/io) are given by perturbative series in a{fio) containing no 
large logarithms. They contain all the information about physics at the scale M which 
is important for low-energy processes. The Wilson coefficients Ci(/i) at low normalization 
scales /i are obtained by solving the RG equations (13.281) with the initial conditions at 
/io ~ M given by matching. The effective theory knows nothing about M; the only 
information about it is contained in Ci(/i). When the effective Lagrangian is applied to 
some physical process with small momenta pi <^ M, it is most convenient to use /i of the 
order of the characteristic momenta: then the results will contain no large logarithms. This 
solution of the RG equation sums large logarithmic terms in perturbation series. 
As already discussed, the anomalous dimension matrix has the block structure 
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in terms of physical and evanescent operators. This means that the evolution equations 
for Wilson coefficients of physical operators don't involve those for evanescent operators. 
The later Wilson coefficients are non-zero; but this does not matter, because renormalized 
evanescent operators vanish. 

As we have seen in Sect. 12.31 the term in the muon loop fl2.28p leads to the contri- 
bution 



AL = cO, c=-^^^ + 0{a^), = {d^F,^){d,F'n (3.29) 



2 a 1 

to the Lagrangian. Due to the equation of motion 



O = e^Oi , (3.30) 
where Oi is defined by (13.21) . Therefore, 

2 a\M) + 0{a^) 

ci(M) = -- — . (3.31) 

We can also understand this more directly by matching the on-shell electron-electron 
scattering amplitudes: 





(3.32) 



g2 3 (47r)'^/2 ^ ' 10 15 {AnY/^ M^+^' ' 

In order to obtain the structure 7(3) ® 7(3), we should have at least three 7 matrices along 
each electron line. Such diagrams first appear at two loops (Fig. [TTl) . Therefore, 

e,(M)^Ml. (3.33) 




Figure 11: A diagram having the 7(3) ® 7(3) structure 
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It is convenient to rewrite the RG equation (13.281) in the form 

dc l'^{c() 



dloga 2/5(a) 



c . 



(3.34) 



At tlie leading (one-loop) order 



The solution of this leading-order equation can be written as a matrix exponent: 

, , /a(/i)\-^o"/(^*) ,,,, 
= -7T7^ c(M) . 



a(M) J 
7oVi = XiVi , 



(3.35) 



If eigenvectors Vi of 7^, 

form a full basi^, the solution (I3.35P has the form 

-M = E-4.(^j 

where the constants Ai are determined by the initial conditions: 

c(M) = J2 ^i^i ■ 

In our particular case 
and the solution of the leading-order RG equation with the initial condition (I3.3ip is 



(3.36) 



(3.37) 



C3(/i) 



1 a\M) 
15x/4l M2 



«M y 

a{M) J 
a{M) 



(2+^41) 



2 + V41 
1 



2 - V41 
1 



(3.38) 



(recall po = -4/3). 



^In some rare cases, there are not enough eigenvectors — Jordan blocks of sizes > 2 appear. It is not 
difficuh to solve RG equations in such cases, too. 
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3.4 Decoupling: photon field 

Now we shall discuss the relation between the full theory and the effective one more sys- 
tematically. Our full theory is QED with massless electrons and muons having mass M. 
When we consider processes with characteristic energies E <^ M, the existence of muons 
is not important. Everything can be described by an effective low-energy theory, in which 
there are no muons. In other words, muons only exist in loops of size ~ if we are 

interested in processes having characteristic distances much larger than 1/M, such loops 
can be replaced by local interactions of electrons and photons. 

The effective low-energy theory contains the light fields — electrons and photons. The 
Lagrangian of this theory, describing interactions of these fields at low energies, contains 
all operators constructed from these fields which are allowed by the symmetries. Operators 
with dimensionalities di > 4 are multiplied by coefficients having negative dimensionali- 
ties; these coefficients contain 1/M*~^. Therefore, this Lagrangian can be viewed as an 
expansion in 1/M. The coefficients in this Lagrangian are fixed by matching — equating 
5*- matrix elements up to some power of Pi/M. 

Operators of the full theory are also expansions in 1/M, in terms of all operators of the 
effective theory with appropriate quantum numbers. In particular, the bare electron and 
the photon fields of the full theory are, up to corrections, 

^o-{ClY^"^'o, Ao={aY^'A', (3.39) 

(primed quantities are those in the effective theory). The bare parameters in the La- 
grangians of the two theories are related by 

eo={Cf'e',, ao = C>o. (3-40) 
The MS renormalized fields and parameters are related by 

where 



(3.41) 



Z,(a(;.), «(;.)) ^^^^^ ' U^^'^ ^^^^^ - W^^"' ^ 



The photon propagators in the two theories are related by 



D±{p){g^.u-^j+aoj^ 



- r° 



-ml t 2\ I PfiPv \ I PfiPu 



(3.43) 



This explains why the same decoupling constant C,a describes decoupling for both the 
photon field A and the gauge-fixing parameter a. It is most convenient to do matching at 
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— > 0, then the power-suppressed terms in f l3.43p play no role. The full-theory propagator 
near the mass shell is 

yos 1 

DAp')-f, Z? = r^. (3.44) 

Only diagrams with muon loops contribute to n(0), all the other diagrams contain no scale. 
In the effective theory 

Z,l(/) = _^. zr = ^-^ = l. (3,45) 

because all diagrams for 11' (0) vanish. Therefore, 

The calculation is essentially the same as in Sect. 12.41 At one loop the bare decoupling 
coefficient is given by (12.531) . For generality, let's suppose there are ni light lepton "flavours" 
in the effective theory and Uf = rii + 1 "flavours" in the full one. Then 

,3.47) 

because n(j9^) is proportional to the number of flavours. With this accuracy, we may put 
a'(/i) = a(yu), and the renormalized decoupling coefficient is given by (I2.56p . It does not 
depend on rii (the case n; = has been considered in Sect. 12. 4p . 

In order to calculate (a (13.421) with two- loop accuracy we need Za and Z^. What is the 
MS renormalization constant Za with nf flavours? It does not depend on masses. If we 
suppose that all flavours have the same mass M, then 11(0) is given by the formula similar 
to (I2.57P where both the one- loop term and the two-loop one are multiplied by ri/ (starting 
from 3 loops, terms with different powers of appear). Re-expressing via renormalized 
quantities at = /iq (I2.67P (with Za = Z^ , see (13.470 . and Z^ (I2.60p ) we have 

(z?)-^.-n(0)^i.l„,^^.(f„,..)(f^)' 

(unneeded higher powers of e are omitted). The inverse quantity Z°l should be equal to 
the minimal renormalization constant ^^(^(/io)) times an expression finite at £ — > 0, and 
we obtain 

4 a; / a \ 2 , ^ 

ZA{a) = l--nf—-2enA — ) . (3.48) 



3 ' Aire ' KAire^ 

The renormalization constant in the effective theory contains ni instead of and 
a' instead of a, where a(/i) and a'{fi) are related by (I3.4ip . and CaifJ') = Cj^if^) (this 
statement will be established in Sect. 13. 6t at the moment, we only need it at one loop, 
see (ESnD). 

With two- loop accuracy, the bare decoupling coefficient fl3.46p is given by (12.580 . be- 
cause only diagrams with a muon loop contribute to 11(0). The renormalized decoupling 
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coefficient fl3.42p is given by fl2.62p : it does not depend on rii (the case = has been 
considered in Sect. 12.41) . Starting from three loops, n;-dependent terms appear. 
The renormahzation group equation 

^^^f^ + 7a(«(^)) - l'Aa\^^)) = (3.49) 
a log /i 

can be used to find /^-dependence of the matching coefficient. As discussed in Sect. 12. 4^ it 
is better to impose the initial condition at /i = M, or at any fi which differs from M by 
corrections of order a. Then (a at this initial /i is 1 + 0{a'^), i.e. corrections are small. 
One popular choice is /iq (12.671) . with Ca(/^o) given by (12.681) . 

3.5 Decoupling: electron field 

The electron propagators in the full theory and in the low-energy theory are related by 

^S{p) = Cl^S'{p) + o(^-^'^ . (3.50) 

It is most convenient to do matching at p — > 0, where power corrections play no role. The 
full-theory propagator near the mass shell is 

1 

where the self-energy of the massless electron is S(p) = Sy(p^)|^. Only diagrams with 
muon loops contribute to Sy(0), all the other diagrams contain no scale; such diagrams 
first appear at two loops (Fig. [T^ . In the effective theory 

SW = f, zr = ^-^ = l. (3.52) 

because all diagrams for Sy(0) vanish. Therefore, 

Z°^ 1 

= = i-Sy(o) ■ ^^-^^^ 




Figure 12: Two- loop diagram for Sy(0) 
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At two loops (Fig. [12]) 

^du If. \ rA / „-\ 2 



- ^^^v{p') = J (l^^^oT'^^ (l^^eoT'^ (^) '(^'9,u - k,KMk') , (3-54) 

where i{k'^g^y — k^ky)Ii{k'^) is the muon-loop contribution to the photon self-energy (12.261) . 
It is transverse; therefore, longitudinal parts of the photon propagators (~ kakp) do not 
contribute, and the result is gauge invariant. We only need the linear term in p in both 
sides: 



i>i:y{Q) = -lel j ^j>^(e^-2p.k^)Yieg,,-k,k^'^^^'^ 



{k 



2\4: 



Averaging over k directions hj p ■ k ^ (k'^/d)^, we obtain 

2(d-l)(rf-4) f cfik n(A;2^ 



Let's calculate the integral 



A. n(*^) _^e|Mqp^__^ (3,56) 



{2tiY {-k'^Y (4vr)"' 
Using (12.27P (and setting Mq = 1) we can reduce it to the vacuum integrals (12.131) (Fig. [5]): 

2 



" d-1 

In particular, 



[AVil, l,n + l)-{d- 2)V{1, 1, n)] . (3.57) 



^' \d-2)id-5Yd~7) - ^^-''^ 



Therefore, we obtain from (13.531) 



etM^'^ 2id - lYd - AYd ~ Q) 
(AttY ^ 'd{d-2){d~5){d-7) 



The renormalized decoupling coefficient is 

cM = ci ^\ (3.60) 

Z^(a(/i),a(/i)) 

Its /i-dependence can always be found by solving the RG equation. It is sufficient to obtain 
it at one point, at some specific /i ~ M, to have the initial condition. The most convenient 
point is /i = M, because a(M) = a'(M) + 0{a^) (Sect. [M]) and a(M) = a'(M) + 0{o?) 
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(Sect. I3.4p . and the differences can be neglected witli our accuracy. Tlie renormalization 
constant up to two loops lias the form 

1 1 / \ ^ 

Z^{a, a) = 1 - ifit^o^ + g [7^0(7^0 + 2/?o) + 7ao7^Io« - 27^1^] (4^) + ' ' ' (3-61) 

where 

is the anomalous dimension of the electron field, 7^0 = 7^,0 + 7^o"' ^^"^ 7^ anomalous 
dimension of the photon field. In QED with nj lepton flavours 

7^(a,a) = 2a— -(4n; + 3)(— ) +■■■ (3.62) 

The effective theory renormalization constant Z'^ is given by a similar formula with primed 
coefficients. Their ratio is 

1 + ^ (7^oA/3o + \MAoi;,a - M^,e^ (^)' , (3.63) 



Z;(a',a') 
where 

4 8 
A/3o = -- , Mao = 3 , = -4 

are the single-flavour contributions to Po, ■jao, li^i- We obtain 

Z;(a,a) ^^UTTe) ' 
Re-expressing fl3.59p via the renormalized a{M), 



we finally obtain 



The RG equation 



dlogfi 



+ 7v(«(/^), «(/^)) - 7^(a'(/^), a'(/^)) = (3.65) 



can be used to find Ci^if^) ^oi 7^ M. In contrast to the case of Ca(/^) fl3.49p . now 7^ — 7^ 
is of order a^, so that changes of of order a (such as, e.g., (12.691) ) don't change the 
coefficient of in (13.640 . 
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3.6 Decoupling: electron charge 



The proper vertex eoT with the external propagators attached (two electron propagators 
5* and one photon propagator D) is the Green function of the fields ipo, ipo, Aq (i.e., the 
Fourier transform of the vacuum average of the T-product of these three fields). Therefore, 
the relation between this quantity in the full theory and in the low-energy effective theory 
is 

eoTSSD = Cj (Ca) e'.T'S'S'D' , (3.66) 
or, taking into account S = C'^S', D = Ca-D', 

eoT'' = iC)-'^' e^" (3.67) 

The vertex r'^(p, p') = 7^ + A^{p,p') on the mass shell has two 7-matrix structures when 
sandwiched between physical spinors; when q = p' — p = 0, only one of them remains: 

= Z^'-ff" . (3.68) 

Only diagrams with muon loops contribute to A^{p,p), all the other diagrams contain no 
scale; such diagrams first appear at two loops (Fig. [T3|) . In the effective theory Z^^ = 1, 
because all diagrams for A'^(p, p) vanish. Therefore, 

Y^ = (^Ot"^^ C{^ = -£^ = Z-, (3.69) 

and we obtain from (13.671) 

c = {cKr'{ar' ■ (3.70) 




Figure 13: Two-loop diagram for A(0, 0) 
The situation in QED is simpler, due to the Ward identity 



Near the mass shell = 0, 



r>.p) = ^. (3,71) 

dp,. 



S{p) = ^, (3.72) 
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and therefore 

Z^'Z^' = 1 . (3.73) 

Similarly, in the effective theory, Z^'^Z'^^ = 1 (in fact, these two renormalization constants 
are equal to 1 separately). Therefore, 

cw = i, c = {cr'. (3.74) 

The bare propagator and vertex are related to the MS renormalized ones by S{p) = 
Z^Sr{p), = Z-pT'f, where Z^, Zp are minimal renormalization constants, and 5*^, are 
finite at £ — > 0. The Ward identity fl3.7ip implies that Z^Z^ is finite at £ — > 0, but the 
only minimal renormalization constant with this property is 1: 

ZrZ^ = 1 . (3.75) 

Therefore, Za = {ZyZ^)""^ Z^^ = Z~^\ similarly, Z'^ = Z'^^, and we obtain 

a = C°f =(C^§)" = C/. (3.76) 

This means that the running charge in full QED (with both electrons and muons) at 
fi = fiQ (12.671) is slightly larger than in the low-energy effective QED (with only electrons): 

«(Ato) = Ca(/^o)a'(/^o) , Ca(/^o) = ^ + y (~^^) ^^''^'^^ 

(see f l2.68p ): the same is true at = Mqs (E 



3.7 Decoupling: bilinear electron currents 

Various operators in the full theory can be expressed via operators in the low energy ef- 
fective theory. All operators with appropriate quantum numbers appear; contributions of 
higher-dimensional operators are suppressed by powers of 1/M. This 1/M expansion of a 
full-theory operator means that its on-shell matrix elements with light external particles 
having small momenta Pi (and physical polarizations), being expanded in pi/M up to some 
finite order, coincides with the 1/M series in matrix elements calculated in the effective 
theory. As an example, here we shall discuss fiavour-nonsinglet bilinear light-fermion cur- 
rents (lA.ip . The full-theory renormalized current is related to the corresponding effective- 
theory current by 

jM = CMjM + Oil/M') . (3.78) 

The on-shell matrix element of jnifJ') is Mn{p,p'; /i) = Z!^^ Z~l {a{fi))Tn{p, p') ■ It should 
be equal to Cj„(/i)M^(p,p'; /"), where M^{p,p';^) = Z'^''Zj~^{a'{fi))T'^{p,p'). Both matrix 
elements are UV-finite; their IR divergences coincide, because both theories are identical 
in the IR region. Any on-shell momenta p, p' can be used; it is easiest to set p = = 0, 
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thus excluding power-suppressed terms. The proper vertex Tn{p,p') = 7(n) + ^n{p,p') at 
p = p' = has the structure r„(0,0) = 7(n)r„, where 7(„) is the Dirac matrix (]A.2p . and 
r„ = 1 + A„ is scalar. Therefore, 

Zj^{a'{fi)) OS 

CjnifJ') = 7 / NN (jn ) Cjn = ^th'^n , (3.79) 

because in the effective theory Z'^^ = 1, = 1. This can also be understood by comparing 
Green functions of ipo, ipo, and jno, similarly to Sect. 13.61 

Only diagrams with muon loops contribute to A„(0,0); the only relevant two-loop 
diagram is shown in Fig. [TH It is equal to 



A„(0,0) = -tel j 



{27tY {k 



2^4 



where Il{k^) is the muon loop contribution to the photon self-energy. Averaging over the 
directions of k and using (]A.6P we obtain 

2 f{d-2nf \ f d^k n(A;2 
A„. = -iel - 1 ' ' 



d J J {27iY{-k^y ■ 

Using (13.591) and the integral (13.581) we obtain 



eX"'V 2. ^ 8(d-6)(r^-l)(r^-rf + l) 

- ^ " 74^^ did-2)id-5)id-7) ■ ^^-^^^ 



The ratio Zjn/Zj^ is given by the formula similar to (13.6^ 

1^ = 1 + i(n - 1) [6{n - 3) - (n - 15)e] (£-) 

and we arrive at [9] 

2 



QniM) = 1 + l(n - l)(85n - 267) (^^^ + • • • (3.81) 




Figure 14: Two-loop on-shell matrix element of a QCD bilinear quark current 
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The three-loop correction has been calculated in [TO] . 

The vector current has Qi = 1 to all orders. For the vector current with a diagonal 
flavour matrix r, the integral of its 0th component is an integer — the difference between 
the numbers of light leptons and antileptons weighted by the diagonal elements of r. This 
difference is the same in the full QED and in the low-energy effective theory. The same 
holds for non-diagonal r by flavour symmetry. We can also see this explicitly. Multiplying 
the Ward identity 

(9S(p) 



r^(o,o)=7 



dp. 



= 7^(1-Sy(0)) 
p=0 



by = [1 ~ ^v(0)]~ , we obtain just Fi = 1. Taking account of the fact that the vector 
current does not renormalize [Zji = 1, Z'^^ = 1) yields Cjiif^) = 1- 

The currents j4 and ja differ from jo and ji by insertion of 7^^. They are related to 
those containing 7^^ by (lA.lSp and (IA.19p . Inserting 7^^ does not change the decoupling 
coefficient. Therefore, 

Z' Z' 

O4 = Cjo^ , Cj3 = Cji-^ , (3.82) 

where Zp^A are given by ( 1A.22I) at two loops, and contain rii instead of n/ = rii + 1. 

The correction to the full-theory vector current ji (see (13.781) ) has been calculated 

in [Tl]. Only one dimension-5 effective-theory operator appears in this correction, namely, 
d'^j[. The coefficient of this operator contains log(yU^/M^). The leading effective-theory 
current j[ also depends on /i when we take into account the leading term ciOi in Li (13.251) : 
its derivative in log/x is proportional to d'^j[. Of course, these two sources of /i-dependence 
cancel, and the full-theory vector current ji does not depend on fi also at the level. 



3.8 Decoupling: electron mass 

In the previous Sections, we considered QED with massless electrons and heavy muons 
(with mass M). Now let's take the electron mass into account as a small correction. We 
shall expand everything up to linear terms in m. The electron propagator in the full theory 
is given by (1A.14P : in the low-energy effective theory, it involves Sy^ instead of T^y^s- These 
two propagators are related by 

1 1 



Comparing the overall factors, we recover (I3.53p . The denominators should also coincide: 
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This means that the on-shell masses are the same in both theories (mos = ''^os) — they can 
be measured at large distances, by a macroscopic experiment (weighting electrons)!^: 

C=!!!^=(C;)-'i±|i(P^. (3.85) 

rriQ ^ ^' 1 + i^skP ) 

This equation should hold for all m <^ M, p <^ M. It is easiest to set m = 0, and 
use (lA.lSp . Then, setting p = 0, we obtain 

C={C-,Y'- (3.86) 

Recalling ( ]A.16p . we finally arrive at 

CM = ^C = go\^)- (3.87) 

In other words, m{^)[^Til)\^ = m' {fi')[iljTilj]'^, does not vary with fi and fi', and does not 
change when one goes from the full theory to the low-energy one. The MS mass decoupling 
constant is 

CmiM) = l- — {^-^] +■■■ (3.88) 



18 V 47r ^ 

its /i dependence can be found from the RG equation 

dlog/i 

(the difference 7^ — 7^ is of order a^). 
3.9 Electron magnetic moment 

When the electron mass m 7^ 0, the electron helicity is no longer conserved. Therefore, the 
effective Lagrangian contains the contribution 

1 

4 

of the magnetic operator Om having dimension 6. 

Let's discuss scattering of an on-shell electron (with a physical polarization) in an 
external magnetic field in full QED (we'll return to the low-energy effective theory later). 
It is described by two form factors ^1^2(0'^). The anomalous magnetic moment is given by 
= -^2(0). In order to find it, we need to expand the vertex function in q up to the linear 
term: 



Lm = --CmOm , Om = metp F^^a'"' ip (3.90) 




tea 



r^ + rf^ + ... . (3.91) 

m 



mv 



^This can also be written as Cm = '^ml'^'m where toq = Z,°^mos, ra'^ = Z'°^ra'^^. 



34 



Then (see, e. g., [T]) 



d-2 



2 1 



(3.92) 



As expected, at the tree level (Fq = 7^, F^"^ = 0) the anomalous magnetic moment is 
/i = 0. 

The first non- vanishing contribution appears at one loop (Fig. [T^ : 



where dots mean higher-loop contributions. 



47r 



(3.93) 




Figure 15: The one-loop contribution to the electron magnetic moment 

There are contributions with muon loops. They first appear at two loops (Fig. [T6ll . 
We shall consider this diagram using the method of regions, see the textbook [1]. In fact, 
in simple situations (like the large-mass expansion) this method is known for a long time, 
see, e.g., |I2]. 




Figure 16: The muon-loop contribution to the electron magnetic moment 
Using the projector (13.921) we obtain 

d'^k 



ten 



d-1 



+ 



(27r 



1 f{d+l){d-2) d'^-d-3 {d + 2){d-3)D2 Dl 



2Do 



+ 



d-2 



d^-4:d + 5 2d^-9d + 13 2{d - 3)D2 



D1D2 



2Dl 
{d-2)Dl 



+ 



(3.94) 
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where the denominators are 



Di = rn? — {k + mv)' 



Do 



and n(A;^) is the muon-loop contribution to the photon self-energy. 
There are two integration regions in the diagram Fig. 



• The muon loop is hard (momenta ~ M), the other loop is soft (momenta ~ m); 

• both loops are hard. 

In the first region, we can expand the muon loop Il(k^) in the small momentum k. By 



dimensionality, these one-loop terms with a sing le scale Mq have the form elM^^'^k'^/M^)'^, 
see (12.281) . When we substitute this expansion into the soft loop integral (with a single 



scale m), these terms produce egM^ 
Fig. M (I333])0, we have 



-2e 



m 



-2ei 



m? /Mq)^. Adding this contribution to that of 



e^ra 



(47r)'^/' 



ci-3 



1 - - 



-2e 



3 (47r)'^/2 



Tie) 1 



(rf- 3)(rf-4) 
5(rf-5) mJ 



+ 



(3.95) 

where dots mean all the other contributions. 

In the purely hard region, the electron lines carry large momenta ~ M, and we can 
expand their propagators in m using Di = D2 — 2mv ■ k. This is a regular expansion (only 
integer powers of appear); by dimensionality, the two- loop integrals (having a single 
scale Mo) produce e^MQ^%m^/M^)''. Avera ging over k directions we see that the first 
non- vanishing contribution appears at 0{m?): 



■Mi 



{AtiY d{d + 2) 



-32 



-4e 



Mi 



32 / a \ 2 
'45M2 [4^ 



(rf-4)2(d-8) 



d{d + 2){d- 2){d -5){d- 7){d - 9) 



(3.96) 



where I3 is defined in (13.561) . 

How can we understand this in the framework of the low-energy effective theory? The 
one-loop contribution (Fig. [15]) is given by the formula (13.931) with Cg — > Cq. The leading 
correction in (13.951) transforms the full-QED charge Cg in the leading term into the effective- 
theory one Cq = {Ca)~^^Q- The contribution of the O^k"^) term in Il{k^) in (13.951) . 



4 e^m"2^M, 



15 M2 



(4vr) 



2e 

° T\e){d-4y 



16 m 



15 M2 V47r/ 



a y 



(3.97) 



-•^^It can be reproduced from p.94p if we formally put n(A;^) = 1, because the photon propagator with 
this insertion becomes simply the Landau-gauge free propagator. 
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corresponds to the contact-interaction diagram (the second diagram in Fig. [TOj) in the 
effective theory (the first diagram in Fig. [10] does not contribute to the anomalous magnetic 
moment, because it has the pure j'^ Dirac structure). Finally, the hard contribution (13.96^ 
is produced by the local vertex of the operator (13.901) in the effective theory: fih = Cm'm? ■ 
Now we see why iih starts from 0{rn?). 

This correspondence between contributions of various integration regions in diagrams of 
the full theory (which contains several widely separated energy scales) and the low-energy 
effective theory (which contains higher-dimensional local operators in its Lagrangian) is 
quite general. Let's briefly discuss one more example — electron-electron scattering at 
small external momenta pi ~ m. In the full theory there are diagrams with muon loop(s). 
For example, let's consider the diagram of Fig. [Ill It contains two integration regions: (1) 
the muon loop is hard, the other one is soft; (2) both loops are hard. In the first region, we 
expand the muon loop in k"^ (cqM^'^^ {k"^ /Mq)"-); the remaining soft loop is a complicated 
non-analytical function of the soft variables: m and Pi. The 0{1) term in Il{k^) just 
replaces Cq — >■ e'^ in the lower-order diagram (Fig. [T7b). The 0(A;^) term corresponds to 
the effective-theory diagram with the 0{a'^) term in the contact interaction ci (Fig. [T7b). 
In the purely hard contribution, the integrand can be expanded in the small variables m 
and Pi. Truncating this regular expansion, we obtain a polynomial in pi which gives a local 
vertex (produced by a local operator in the effective Lagrangian), Fig. [TTb . In particular, 
the zeroth term of this expansion produces 0{a^) contributions to the Wilson coefficients 
ci and C3. 




a be 



Figure 17: Contributions to electron-electron scattering 



4 Decoupling of heavy flavours in QCD 
4.1 Effective low-energy QCD 

In QED, effects of decoupling of muon loops are tiny. Also, pion pairs become important 
at about the same energies as muon pairs, so that QED with electrons and muons is a 
model with a narrow region of applicability. Therefore, everything we discussed in Sect. [3] 
is not particularly important, from the practical point of view. 

In QCD, decoupling of heavy flavours is fundamental and omnipresent. It would be 
a huge mistake to use the full 6-fiavour QCD at characteristic energies of a few GeV, or 
a few tens of GeV: running of as(/i) and other quantities would be grossly inadequate, 
convergence of perturbative series would be awful because of large logarithms. In most 
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cases, anybody working in QCD uses an effective low-energy QCD, where a few heaviest 
flavours have been removed. Therefore, it is important to understand decoupling in QCD. 
And to this end the lessons of QED are very helpful. 

Let's consider QCD with a single heavy flavour having mass M; for simplicity, all other 
flavours are supposed to be massless. Then the behaviour of light quarks and gluons at 
low momenta -C M is described by the low-energy effective theory. Its Lagrangian is the 
usual QCD Lagrangian (of course, without the heavy-quark field) plus higher- dimensional 
terms (whose coefficients are suppressed by powers of 1/M). Power corrections to the 
Lagrangian first appear at dimension 6. 

Unlike the QED case (12. 2p . there exists a non- vanishing dimension-6 operator con- 
structed from 3 G^j^: 

Ol, = 9or''Gl/Gl;Gl,^ (4.1) 

(a similar operator with d"'^'^ is identically zero, as in QED (12. 2p ). One- loop renormalization 
of the operator (14. ip has been considered in [131 Ej- There is also the operator 



0% = {D^^G^,,^)iD,Gfn (4.2) 

which reduces to a quark operator due to equations of motion. The operator G^^^D'^Gq^^ 
reduces to (14. 2p . similarly to the QED case (Sect. 12.11) . 
The local four-quark dimension-6 operators are 



^ q ' ^ q ' 

= (Y ^07{n)t"gO^ (Y ^07(n)t"gO^ 



(4.3) 



Only operators with odd n conserve the light-quark helicity. Those with even n can only 
appear in the Lagrangian being multiplied by a power of light-quark masses; these terms 
have higher dimensions, and we don't consider them here. Operators (14. 3 p with n > 5 
are evanescent, and are eliminated by renormalization (Sect. 13. 2p . Only 4 operators Ogi, 
Oq3, Oqi, Ogs appear in the effective QCD Lagrangian (at the dimension-6 level). The 
operator — vanishes due to equations of motion. Renormalization of the four- 

quark operators (14. 3 p can be considered similarly to Sect. 13. 2t the only extra complication 
is their colour structures. 

The 0{q'^) term in the one-loop gluon self-energy (I2.26P (which is given by (12.280 with 
Cq — > Qq and the colour factor Tp) gives the term c°20°2 in the Lagrangian with 

.»,(M)--^f (^ + 0(a.;(M))) (4.4) 

(of course, our one-loop calculation leaves higher orders in unknown; if we choose /z ~ M, 
there are no large logarithms, and the error due to truncation of the perturbative series is 
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minimal). We are interested only in the S'-matrix; therefore, we can eliminate this term in 
favour of Cg^O^^ with 

5.1 (M) = {al{M) + 0{al{M))) . (4.5) 

The 3-gluon interaction c°]^0°^ also first appears at one loop. The coefficient can 
be found by matching the 3-gluon scattering amplitudes (Fig. [18], pi + P2 + = 0) in the 
full theory and in the effective one. Diagrams without heavy-quark loops are the same in 
both theories. We are interested only in diagrams with (at least one) heavy-quark loop in 
the full QCD. Contributions of the fully hard regions in such diagrams can be expanded 
in the small momenta pi up to cubic terms. We cannot put all 3 legs on-shell (as we did 
in the case of the 4-photon scattering amplitude in Sect. 12. 2p . Therefore, when doing the 
effective-theory calculation, we have to take both operators (14. ip . (14. 2 p into account. 




Figure 18: Three-gluon scattering amplitude 



The scattering amplitude appears at the tree level — just the 3-gluon vertex. It should 
be multiplied by the one-loop external-legs renormalization factor Zj( . All contributions 
to Za without a heavy-quark loop appear also in the effective theory and hence cancel. 
Leaving only the heavy-quark contribution to Za we get 



(4.6) 



where the 3-gluon vertex tensor structure ^^1^*2^3 jg linear in the external momenta p,. 
There are also one-loop vertex diagrams with a heavy-quark loop: the diagram in Fig. [19] 
and the one with the opposite direction of the quark line. The d°-^°-^°-^ contributions cancel 
(just like in QED) and the /"i"2a3 ones are equal (so that we can calculate the f"-^"-'^"-^ term 
of the diagram in Fig. [19] and then double it). We expand this diagram in the external 
momenta pi. Only odd powers appear; we need terms up to cubic ones. Averaging the 
integrand over k directions and using (I2.10p we obtain 



-2e 



0-10-20-3 



(47r)^/2 



4 d — A 

yMl/^2M3 + ^ 

3 180 



(4.7) 



where the 3-gluon matrix elements f^-f^^asj^t^u^ifj-s ^^iq operators are cubic in pi. 
Adding (14.61) and (14.71) we see that terms linear in pi produce —g'^f^-^^-^^-^v^^^f^^^^^, the 
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Figure 19: Heavy-quark loop contribution to the 3-gluon scattering amplitude 

elementary 3-gluon vertex in the effective theory. The structure T^^^^''^ gives c^^, the 
coefficient of the 3-gluon operator (14. ip in the effective Lagrangian [15] 



=»'(-^)-o^(^-^(«H- ''''' 

The structure T^'i'^^/is reproduces (14. 4p . It is not difficult to extend this analysis to two 
loops using (I2.13p . 

4.2 Decoupling 

Running of the full-theory coupling Q;i"'^^^(yu) is governed by the (ni + l)-flavour /3-function; 
running of the effective-theory coupling ai"'''(/i) is governed by the n^-flavour /3-function; 
their matching is given by 

= C«(/^)«i"'H/i) , (4.9) 

with [IS] 

C.(/^o) = 1 + f^C. - ^-C^ T, { + . . . (4.10) 



.3 9 y V 47r , 

where /xq is defined by f l2.67p . Here the Cp term can be obtained from the QED result (I3.77P 
by inserting the obvious colour factors; the Ca term is more difficult to obtain. The three- 
loop result has been obtained in [3] , and the four- loop one in [TTJ [Tsj"] . The RG equation 

"^^"f^"^^^ + 2/3("'+i)(a^+i)(/i)) - 2/5("')(a^)(/i)) = (4.11) 
ct log /i 

can be used to find Ca(/i) for /i 7^ M. 

The QCD running coupling a^i^^i) not only runs when /i varies; it also jumps when 
crossing heavy-flavour thresholds. The behaviour of as(/i) near is shown in Fig. [20] (this 
figure has been obtained using the Mathematica package RunDec [20], which takes into 
account four-loop /5-functions and three-loop decoupling). At /i > mb, the correct theory is 



^^The result of |18| contains one master integral which was not known analytically, only numerically, 
with 37-digits precision. An analytical expression for this integral has been published later [19] . 
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0.225 n 




0.21 -I > , 

nib — O.SGeV m^, + 0.5GeV 



Figure 20: af\fi) and aT'in) 



(4)/ 



the full 5-fiavour QCD {a^^\fi), the solid line); at /i < rrib, the correct theory is the effective 
low-energy 4-fiavour QCD {ai^\ij,), the solid line); the jump at fi = rub f l4.10p is shown. Of 
course, both curves can be continued across (dashed lines), and it is inessential at which 
particular n ~ rUb we switch from one theory to the other one. However, the on-shell mass 
(or any other mass which differs from it by 0{as), such as, e.g., /xq) is most convenient, 
because the jump is small, C(af). For, say, = 2mb or /i = mb/2 it would be O^a"^). 

Light-quark masses mi{fi) also rum with fi, and also jump when crossing a heavy-quark 
threshold. The QCD result 



m 



(M) =m('^')(M) 



89 



CpTp 



as{M) 
47r 



+ 



(4.12) 



can be obtained from the QED one (13.881) by inserting the obvious colour factors. 



5 Conclusion 

In the past, only renormalizable theories were considered well-defined: they contain a finite 
number of parameters, which can be extracted from a finite number of experimental results 
and used to predict an infinite number of other potential measurements. Nonrenormalizable 
theories were rejected because their renormalization at all orders in nonrenormalizable 
interactions involve infinitely many parameters, so that such a theory has no predictive 
power. This principle is absolutely correct, if we are impudent enough to pretend that our 
theory describes the Nature up to arbitrarily high energies (or arbitrarily small distances). 
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Our current point of view is more modest. We accept the fact that our theories only 
describe the Nature at sufficiently low energies (or sufficiently large distances). They are 
effective low-energy theories. Such theories contain all operators (allowed by the relevant 
symmetries) in their Lagrangians. They are necessarily nonrenormalizable. This does not 
prevent us from obtaining definite predictions at any fixed order in the expansion in E /M, 
where E is the characteristic energy and M is the scale of new physics. Only if we are lucky 
and M is many orders of magnitude larger than the energies we are interested in, we can 
neglect higher-dimensional operators in the Lagrangian and work with a renormalizable 
theory. 

Practically all physicists beheve that the Standard Model is also a low-energy effective 
theory. But we don't know what is a more fundamental theory whose low-energy approx- 
imation is the Standard Model. Maybe, it is some supersymmetric theory (with broken 
supersymmetry); maybe, it is not a field theory, but a theory of extended objects (super- 
strings, branes); maybe, this more fundamental theory lives in a higher- dimensional space; 
or maybe it is something we cannot imagine at present. The future will tell. 

Acknowledgements. I am grateful to K. G. Chetyrkin and M. Misiak for useful 
discussions. 



A Renormalization of bilinear currents 

Here we shall consider the currents 

jnO = V'07(n)TV'0 , (A.l) 

where 

l[n)=l^''---l'-\ (A.2) 

in QED with Uf > 2 lepton flavours, where r is a flavour matrix with Trr = (flavour- 
singlet currents have some peculiarities, because the fermion line emerging from the oper- 
ator vertex can return to the same vertex again; we shall not discuss such currents). The 
renormalized currents are related to the bare ones by 

jnO = Zjn{a{fx))jn{l^) , (A. 3) 

where the Zjn are minimal renormalization constants. The dependence of jn(A*) is deter- 
mined by the renormalization-group equation 

^ +7in(a(/^)) ) Jn(/x) = 0, (A.4) 



dlog/i 



where 

7> = ^ (A.5) 
a log// 

is the anomalous dimension. 
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Let the sum of one-particle- irreducible bare diagrams with a current vertex 7(„), an 
incoming fermion with momentum p, and an outgoing fermion with momentum p', not 
including the external fermion propagators, be the proper vertex r„(p, p') = 7(„) +A„(p,p') 
(Fig. EH). When the vertex is expressed via the renormalized quantities a(/i), it should 
become Zm^niPjP'), where the renormalized vertex r^{p,p') is finite in the limit e — > 0. 
When the proper vertex of the renormalized current Zj'^Tn{p,p') is multiplied by the two 

1 /2 

external- leg renormalization factors , it should give a finite matrix element. Therefore, 
Zjn = ZgZm- The UV divergences of A„(p, p') do not depend on the quark masses and 
the external momenta. Therefore, we may assume that all fermions are massless, and set 
p = p' = 0. An IR cut-off is then necessary in order to avoid IR 1/e terms. 



I I 
I I 




Figure 21: Proper vertex of a bilinear current 
In the one-loop approximation (Fig. [2Tb). 



A„(0,0) = -ze^ 



{27rY (P)3 



d 



where the averaging |^7(n)^ —>■ k'^j'^'^'Ja/ d has been used. Now we use an important identity 

7^7{n)7/.= (-ir(c?-2n)7(„). (A.6) 

It can be understood as follows. Let's take some specific component of 7(„). Then n values 
of the index /i which are present in this component of 7(„) contribute — (— while the 
remaining d — n values of fi contribute {—l)"'{d — n). Using the UV divergence 



d'^k 



{27tY (A;2^2 



UV 



(47r)% 



we obtain 



Zrn = l + £^^[{n-2r-l + a] 



The gauge dependence is cancelled by 



1 — a- 



a 

'4:7ie 



and we arrive at 



a 



7,„ = -2(n-l)(n-3)— + ■ 



(A.7) 
(A.8) 

(A.9) 
(A.IO) 
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Up to two loops [21] 



'jjn= ~2{n-l){n-3)-^ 



l + l(5(n-2f-19)^' 



{A.ll] 



The three-loop result has been derived in |22j . 

The vector current has 7^1 = to all orders. For the vector current with a diagonal 
flavour matrix r, the integral of its 0-th component is an integer — the difference between 
the numbers of light fermions and antifermions weighted by the diagonal elements of r. It 
does not depend on fi. The same holds for non-diagonal r by flavour symmetry. We can 
also see this explicitly. Multiplying the Ward identity 

r^(p,p) = ^^ (A.12) 

by transforms into S~^, and hence makes the left-hand side finite. Therefore, 
Zji = Z^Zyi = 1. 

The renormalization of the scalar current {n = 0) is closely related to that of the 
fermion mass m. The MS renormalized mass m{fi), 

mo = Z,rn{(y{fi))m{fi) (A. 13) 

(with a minimal renormalization constant Zm), is defined in such a way that the bare 
propagator 

" ?i - m„ - E(p) " 1 - Ev(p2) ^ 1 + '^-"^ 

(where the fermion self-energy S(j9) = ^J^vip"^) +^o^s{p'^)), when expressed via the renor- 
malized quantities a(/i), a(/i), m{fi), is equal to Z^Sr{p), where the renormalized propa- 
gator Sr{p) is finite in the limit e — > 0. In order to find Z^, it is sufficient to consider 
Ip^I ^ m^ and retain only terms linear in m. Differentiating diagrams for S in mo, we get 

l+S<j(p2)L^^^ = ro(p,p). (A.15) 

The renormalization constant Z^ is defined by the condition that ZqZm{^ + S5) is finite. 
Therefore, 

Zm = Z-^' . (A.16) 

In other words, m(fi) ['ipT^']^ = fno4'o'T'4'o is not renormalized. The anomalous dimension 
of the mass is 

a 20 \ ( oi\^ 
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It is not possible to introduce an anticommuting 75 satisfying 7^*^7^+7^7^'~^ = in c? di- 
mensions. A way out was proposed by 't Hooft and Veltman. Let us split our ci- dimensional 
space-time into a 4-dimensional subspace and the orthogonal (d — 4)-dimensional subspace, 
and define 7^^^ = — Z7°7"'^7^7'^ = (z/4!)£:q,/3^57°7'^7'''7'^, where the tensor EaP'^s lives in the 
4-dimensional subspace. This 7^^ anticommutes with 7^^ from the 4-dimensional subspace 
and commutes with those from the ((i — 4)-dimensional one. This consistent definition leads 
to the anomaly of the fiavour-singlet axial current. But the fermion line cannot be closed 
in matrix elements of fiavour-nonsinglet currents. Therefore, traces with a single 75, which 
can lead to an anomaly, never appear, and one may use a naively anticommuting 7^*-^ with- 
out encountering contradictions. The pseudoscalar currents Jac(/^) = ZJ'J {a{iJ,))%jjo'y^'~"Tijjo 
and jHv(/^) = Zj'^ {a{fi))ijjo'yf^ TipQ are related to each other by a finite renormalization: 

jAc(/i) = ^p(a(/u))jHv(/i), Zp{a) = 1 + zpi— + zp2 i^—) +■■■ (A.18) 
Similarly, the axial currents are related by 

j^c(/i) = Z^(a(/x))j^^v(/i)- (A.19) 

This is clearly discussed in [23], where references to earlier papers can be found. 

The finite renormalization constants Zp^ can be found by calculating matrix elements 
of the operator equalities flA.lSp . flA.lQD . The easiest choice is the matrix element from an 
initial fermion to a final one with momenta p, p (or p, 0). At one loop (but not beyond) 
Zj4^ = Zjo, Zj3 = Zji = 1; calculating the proper vertices and equating the matrix elements, 
we can easily get 

Zp(«) = l-8^ + ---, Z^{a) = l-A^ + --- (A.20) 
4tt 47r 

They can be also obtained from the anomalous dimensions of the currents. Differentiat- 
ing flA.lSD and flA.lQD . we have 

d\ogZp{a) _ 7jo(a) - 7^4(0) 

rfloga 2/5(a) ' 

dlogZA{a) 7ji(a) -7j3(a) , „ 

dloga 2/3{a) 

The one- loop results ([OOl) can be reproduced using the two-loop anomalous dimen- 
sions (lA.lip . Now we see the reason why the last term in (lA.lip . which is not symmetric 
with respect to n — 4 — n, is proportional to /?o- As Zp A can be obtained solely from the 
anomalous dimensions, they are determined by the UV behaviour of the matrix elements, 
and cannot depend on masses and external momenta. 
The two-loop results 



a o /ax 



2 



Zp(a) = 1 -8— + -Uf { — ] +■••, 

4vr 9 'UJ ' (A.22) 

Z^(a) = l-4- + -(99 + 2n,)(-) +... 
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can be obtained either from two-loop matrix elements or from the three-loop anomalous 
dimensions. The three-loop results have been calculated [23] from the matrix elements 
only. 

Why have we not discussed a similar relation j^Q^ii) = ZT{a{fi))j^{fi) between the 
tensor currents Jq^q = tpo'j^^a'^'^Tipo and Jq^v = "4^01^^ o'^^T'ipa! The current j^^ has 
the same anomalous dimension 7^2 as the current jj*^ without 7^*". Multiplication of a^^ 
by 7^^^ is merely a space-time transformation, e.g., 7|^^cr'^^ = —ia"^^, and hence j^y has 
the same anomalous dimension, too. Therefore, Zxiaifi)) cannot depend on fi, and we 
conclude that Zt^o) = 1. 



References 

[1] A. G. Grozin, Lectures on QED and QCD: Practical Calculation and Renormalization 
of One- and Multi-Loop Feynman Diagrams, World Scientific (2007). 

[2] G. V. Dunne, Heisenherg-Euler Effective Lagrangians: Basics and Extensions, in 
From Fields to Strings: Circumnavigating Theoretical Physics, ed. M. Shifman, 
A. Vainshtein, J. Wheater, World Scientific (2005), v. 1, p. 445 
(8jXiv:hep-th/0406216| . 

[3] K. G. Chetyrkin , B.A. Kniehl, M. Steinhauser, Nucl. Phys. B 510 (1998) 61 
|arXiv:hep-ph7 9708255| . 

[4] V.A. Smirnov, Applied asymptotic expansions in momenta and masses, Springer 
Tracts in Modern Physics 177, Springer (2002). 

[5] H. Euler, B. K5ckel, Naturwissensch. 23 (1935) 246; 

W. Heisenberg, H. Euler, Z. Phys. 98 (1936) 714 | arXiv:physics/ 0605038|; 

V. Weisskopf, Kong. Dans. Vid. Selsk. Math-fys. Medd. XIV No. 6 (1936); in 

A. I. Miller, Early Quantum Electrodynamics: A Source Book, Cambridge University 

Press (1994). 

[6] X.-w. Kong, F. Ravndal, Nucl. Phys. B 526 (1998) 627 [ arXiv:hep-ph/9 803216|. 

[7] M.J. Dugan, B. Grinstein, Phys. Lett. B 256 (1991) 239. 

[8] K. G. C hetyrkin, M. Misiak, M. Munz, Nucl. Phys. B 520 (1998) 279 

|arXiv:hep-ph/9711280 |; Nucl. Phys. B 518 (1998) 473 [arXiv:hep-ph/9711266| . 

[9] A.G. Grozin, Phys. Lett. B 445 (1998) 165 |arX iv:hep-p h/9810358]. 

[10] A.G. Grozin, A. V. Smirnov, V.A. Smirnov, JHEP 11 (2006) 022 
(a^v:hep^h7 0609280| . 

[11] K. G. Chetyrkin, Phys. Lett. B 307 (1993) 169. 



46 



[12 
[13 
[14 

[is: 

[16 

[17; 

[is: 

[19 

[2o: 

[21 

[22 
[23 



S.G. Gorishny, S.A. Larin, Nucl. Phys. B 283 (1987) 452. 
S. Narison, R. Tarrach, Phys. Lett. B 125 (1983) 217. 

A.Yu. Morozov, Sov. J. Nucl. Phys. 40 (1984) 505; Preprints ITEP-190, 191 (1983). 

V. A. Novikov, M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, Fortsch. Phys. 32 
(1984) 585. 

S. A. Larin, T. van Ritbergen, J. A.M. Vermaseren, NucL Phys. B 438 (1995) 278 
(arXiv:hep-ph/9411260|. 

K. G. Chetyrkin, J. H. Kiihn, C. Sturm, NucL Phys. B 744 (2006) 121 
|arXiv:hep-ph/0512060|. 



Y. Schroder, M. Steinhauser, JHEP 01 (2006) 051 [a rXiv:hep-ph/0512058| . 

B.A. Kniehl, A.V. Kotikov, A. L Onishchenko, O.L. Veretin, Phys. Rev. Lett. 97 
(2006) 042001 |arXiv:hep-ph70 607202 1 . 

K. G. Chetyrkin, J . H. Kiihn, M. Steinhauser, Comput. Phys. Commun. 133 (2000) 
43 |arXiv:hep-ph/0 004189| . 

D. J. Broadhurst, A. G. Grozin, Phys. Rev. D 52 (1995) 4082 
|arXiv:hep-ph/9410240j. 



J. A. Gracey, Phys. Lett. B 488 (2000) 175 | arXiv:hep-ph/0007171| . 

S. A. Larin, The renormalization of the axial anomaly in dimensional regularization, 
in Quarks-92, ed. D.Yu. Grigoriev, V. A. Matveev, V. A. Rubakov, P. G. Tinyakov, 
World Scientific (1993) p. 201 [arXiv:hep-ph/9302240|. 



[24] S.A. Larin, J. A. M. Vermaseren, Phys. Lett. B 259 (1991) 345. 



47 



